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The compound KTi(S04)2-H20 was recently reported as a quasi one-dimensional spin 1/2 com- 
pound with competing antiferromagnetic nearest neighbor exchange Ji and next-nearest neighbor 
exchange J2 along the chain with a frustration ratio a = J2/J1 ~ 0.29 [Chem. Mater. 20, 8 (2008)]. 
Here, we report a microscopically based magnetic model for this compound derived from density 
functional electronic structure calculations along with respective tight-binding models. Our calcula- 
tions confirm the quasi one-dimensional nature of the system with antiferromagnetic Ji and J2 , but 
suggest a significantly larger frustration ratio q « 1.1 ± 0.2. Based on transfer matrix renormaliza- 
tion group calculations we found that, due to an intrinsic symmetry of the Ji- J2 model, our larger 
frustration ratio a is also consistent with the previous thermodynamic data. To resolve this issue, 
we propose performing high-field magnetization measurements and low temperature susceptibility 
measurements which should allow to precisely identify the frustration ratio a. 

PACS numbers: look it up 
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I. INTRODUCTION 

For several decades, low dimensional magnetism has 
attracted great interest in solid state physics and chem- 
istry. Since the influence of quantum fluctuations be- 
comes especially pronounced for low dimensional spin- 
1/2 systems, these systems have been investigated exten- 
sively, both theoretically and experimentally. Quantum 
fluctuations become even more important in determining 
the ground state and the nature of the low lying excita- 
tions if the system under consideration exhibits strongly 
frustrated interactions. 

Whereas pure geometrical frustration, i.e triangular, 
Kagome or pyrochlore lattices, can be realized by spe- 
cial symmetries in two or more dimensions, frustration in 
one-dimensional systems (linear chains) is generally real- 
ized by competing interactions. The simplest frustrated 
one dimensional model is the J1-J2 model with nearest 
neighbor (NN) exchange Ji and next nearest neighbor 
(NNN) exchange J2, where J2 is antiferromagnetic. The 
phase diagram of this seemingly simple model is very 
rich. Depending on the frustration ratio a = J2/J1, a 
variety of ground states are observed in corresponding 
quasi-lD systems: (i) ferromagnetically ordered chains in 
Li2Cu02^'2 (0 < |a| < 0.25); (ii) helical order with differ- 
ent pitch angles in LiVCu04, LiCu202, and NaCu202^"~ 
{a < —0.25); (iii) spin-Peierls transition in CuGeO- )^°d^ 
{a > 0.2411), i.e. both exchange couplings are antiferro- 
magnetic. 

The possibility of subtle interplay between spin, or- 
bital, charge and lattice degrees of freedom due to the 
threefold orbital degeneracy of Ti"^"*" in octahedral envi- 
ronments, makes titanium 3+ based oxides an interesting 
class of materials to study. Exotic features like orbital- 
liquid state in LaTiOs and presence of strong orbital 
fluctuations in YTiOa have been reported i^^d^ While an 
abundance of experimental work exists on low dimen- 



sional spin- 1/2 cuprates (with Cu^+), materials based on 
spin-1/2 titanates (with Ti^"*") are rather sparse. A well 
known example of low-dimensional titanates is the new 
class of inorganic spin-Peierls materials TiOX {X = CI, 
Br).^ Quasi one-dimensional magnetism was observed in 
TiOCl and TiOBr along with a first-order transition to a 
dimerized non-magnetic ground state (spin-Peierls like) 
below 67 K and 27 K, respectively. The main difference 
between S=l/2 titanates and vanadates as compared to 
the cuprates (3d^ vs. 3d^) is that the unpaired electron 
resides in the t2g complex for the former, while occupy- 
ing the Cg complex for the latter. This usually results in 
narrow bands at the Fermi level (E^) for the titanates 
and vanadates as compared to the cuprates, which in 
turn leads to small values for the exchange couplings and 
brings several experimental conveniences. Since the tem- 
perature scale for the magnetic contribution to the spe- 
cific heat is of the order of J, separating magnetic and 
phonon contribution in the specific heat (Cp) measure- 
ments is thus relatively easy. For the same reason, mag- 
netization measurements can reach the saturation mo- 
ment in experimentally attainable fields, thereby provid- 
ing more information about the exchange parameters and 
the frustration regime. Thus, in many respects, spin 1/2 
compounds with a singly occupied t2g orbital are an ideal 
object to study the physics of quasi one-dimensional frus- 
trated chains. 

KTi(S04)4-H20, a new member in the family of ti- 
tanium alums has recently been synthesized^^. The ti- 
tanium ions are in the 3+ (d^) oxidation state in this 
material. Specific heat and susceptibility measurements 
indicate that KTi(S04)4-H20 is a S=l/2 frustrated chain 
system. Fits to the susceptibility data using exact diago- 
nalization on up to 18 spins resulted in estimates for lead- 
ing exchanges Ji = 9.46 K and J2 = 2.8 K, both AFM 
with a = 0.29. In the well studied J1-J2 frustrated chain 
model with both interactions being AFM, the model un- 




FIG. 1. (Color Online) Top: Crystal structure of 
KTi(S04)2-H20. Double chains of TiOe octahedra run along 
the b-axis. The octahedra are connected on either sides by 
SO4 tetrahedra. Water molecules are bound to the octahe- 
dra and separating the double chains along the c-axis. The 
potassium atoms (not shown here) separate the double chains 
along the a-axis. Bottom: An isolated segment of the double 
chain. The nearest neighbor (NN) and next nearest neighbor 
(NNN) hopping paths are represented as ti and t2 respec- 
tively. 



dergoes a quantum phase transition at a = 0.2411 to 
a two-fold degenerate gapped phase, and for a > 0.2411 
the system will exhibit spontaneous dimerizationJ^"— In 
light of this result, KTi(S04)4-H20 would be expected 
to occupy the highly interesting region of the phase dia- 
gram, with a small gap A < Ji/20J^i^ 

Here, we report the results of an electronic structure 
analysis from first principles, carried out to obtain a mi- 
croscopic picture of the origin of the low dimensional 
magnetism in KTi(S04)4-H20. The magnetically active 
orbital is identified using band structure calculations fol- 
lowed by subsequent analysis of the exchange couplings. 
The calculated J's have the same sign as the experiments 
{i.e. AFM), but the estimated frustration ratio a is con- 
siderably larger compared to the experimental findings. 
We explore in detail the reasons for this discrepancy and 
propose an alternative solution that fits the experimen- 
tal data, as well as being consistent with our band struc- 
ture calculations. To this end, we have simulated the 
temperature dependence of magnetic susceptibility us- 
ing the transfer-matrix renormalization group (TMRG) 



method to unambiguously identify the exchange cou- 
plings that describe the microscopic magnetic model for 
KTi(S04)2-H20. On the basis of these results, we sug- 
gest performing high-field magnetization measurements, 
which will be a decisive experiment to identify the precise 
frustration ratio a. The influence of crystal water on the 
observed ground state is also discussed in detail. 

The remainder of the manuscript is organized as fol- 
lows: The crystal structure of KTi(S04)2-H20 is de- 
scribed in section II. The details of the various calcu- 
lational methods are collected in section III. The results 
of the density functional theory based calculations, in- 
cluding the band structure and the accordingly derived 
microscopic model is described in section IV. In section 
V we demonstrate the effects of the crystal water in this 
compound by calculating the Wannier functions. The 
outcome of the TMRG simulations is compared to the 
previous experiment in section VI, which is followed by 
a discussion and summary in section VII. 



II. CRYSTAL STRUCTURE 

Throughout our calculations, we have used the re- 
cently determined experimental^" lattice parameters in 
the monoclinic space group (P2i/m) of KTi(S04)2-H20: 
a = 7.649 A, 6 = 5.258 A, c = 9.049 A and ;3 = 101.742°. 
The crystal structure which is isomorphous to that of 
the mineral Krausite^^, KFe(S04)2-H20 is displayed in 
Fig. [TJ Isolated pairs of chains ( "double chains" ) of TiOe 
octahedra run along the crystallographic 6-axis. The oc- 
tahedra are distorted and have no edge sharing or corner 
sharing oxygen atoms. The SO4 tetrahedra corner share 
with three adjacent TiOg octahedra, forming an isosce- 
les triangle. These trangles edge-share to make up the 
double chains. The single chains are displaced with re- 
spect to each other both laterally and vertically. Large 
K+ ions isolate the double chains along a, while along c 
the chains are separated by water molecules that share 
the oxygen atom with the TiOg octahedra. The water 
molecules are oriented in the ac plane. All of the octahe- 
dral 0-Ti-O bond angles deviate slightly away from 90 
and there are three pairs of Ti-0 bond lengths of 2.001 A, 
2.056 A (in the ab plane) and 2.043 A (along c). The 
shortest Ti-Ti distance is 4.93 A, and is between near- 
est neighbors (NN) on the adjacent chains (ii in Fig. [T|). 
Within the chains, neighboring Ti are 5.23 A apart(i2 in 
Fig. [T|). By analogy with the Ji — J2 Heisenberg model, 
we will call the magnetic interactions corresponding to 
the shorter distance (NN) Ji and the longer distance 
(next nearest neighbors, NNN) J2. In case of a perfect 
octahedral environment, the three t2g states would be 
degenerate and thus warrant additional effects {i.e. lat- 
tice distortion, spin/charge/orbital ordering) to lift the 
degeneracy and allow for an S=l/2 singlet ground state 
for a Ti3+ ion,^ In KTi(S04)2-H20, there is a smah 
distortion of the octahedra and consequently splitting of 
the t2g levels can be expected. The related case of Ti- 



OCl, another system containing Ti^+ ion, also possesses 
a distorted arrangement of TiCl2 04 octahedra, though 
the distortions are much larger with equatorial Ti-0 and 
Ti-Cl bond lengths of 2.25 A and 2.32 A, respectively, 
and an apical Ti-0 bond length of 1.95 A. Consequently, 
the t2g orbitals in TiOCl were thought to split into a 
lower energy dxy and higher energy dxz.yz orbitals. Elec- 
tronic structure calculations confirmed this interpreta- 
tion and revealed the magnetically active orbital for the 
S=l/2 chains in TiOCl was indeed the lower energy d^y 
orbitalfii!^^ though a prolonged discussion of possible or- 
bital fluctuations ensued afterwards. Therefore, an anal- 
ysis of KTi(S04)2-H20 from the structural point of view 
alone is not sufficient to determine the ground state of 
the system. Detailed calculations are necessary to un- 
derstand the correct orbital and magnetic ground state 
of the system. 
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III. CALCULATIONAL DETAILS 

The DFT calculations were performed using a full po- 
tential nonorthogonal local orbital code (FPLO) within the 
local (spin) density approximation (L(S)DA)^i21. The 
energies were converged on a dense k mesh with 300 
points for the conventional cell in the irreducible wedge of 
the Brillouin zone. The Perdew and Wang flavor^^of the 
exchange correlation potential was chosen for the scalar 
relativistic calculations. The strong on-site Coulomb re- 
pulsion of the Ti id orbital was taken into account us- 
ing the L(S)DA+[/ method, applying the "atomic limit" 
double counting term. The projector on the correlated 
orbitals was defined such that the trace of the occupa- 
tion number matrices represent the id gross occupation. 
The AFM parts of the exchange couplings arc computed 
as J^^^ = Atj /Ucs, by mapping the results of the LDA 
calculations on to a tight binding model (TBM) which is 
then mapped on to a Hubbard model, and subsequently 
to a Heisenberg model because the system belongs to the 
strong correlation limit C/cff ^ i (t is the leading trans- 
fer integral at half filling). The full exchange couplings 
are obtained by mapping the LSDA-|-C/d total energies 
of various supercells with collinear spin configurations 
to a classical Heisenberg model. The supercells used to 
calculate Ji and J^ were sampled using 300 and 100 k 
points, respectively. Maximally localized Wannier func- 
tions (WF) were calculated for the Ti i2g orbitals, also 
using FPLO^S, to obtain a visual insight of the relevant 
orbitals and superexchange paths. 

The magnetic excitation spectrum of frustrated spin 
chains was simulated using exact diagonalization code 
from the ALPS package. ^° We used periodic boundary 
conditions and considered finite lattices comprising up 
to iV = 32 spins. 

The magnetic susceptibility of infinite frustrated spin 
chains was simulated using the transfer-matrix renormal- 
ization group (TMRG) technique. "^^ For each simulation, 
we kept 120-160 states, the starting inverse temperature 



FIG. 2. (Color Online) Top: Total and partial DOS obtained 
within LDA for KTi(S04)2-H20. The valence panel is pre- 
dominantly comprised of Ti-3d and 0-2p states. The sulphur 
states (not shown here) lie mainly below -8 eV. The contri- 
bution from K, S and H sites are negligible in the displayed 
energy range. Bottom: Ti 3d-orbital resolved DOS. The t2g 
and eg complexes are split by a ligand field splitting of about 
2 eV. The d^z orbital and the dy^ orbital are very close and 
split from the broad dxy (larger bandwidth) orbital. 



was set to 0.05Ji, and the Trotter number was varied 
between 4-10^ and 16-10'^. The results were well con- 
verged for the whole t emp erature range of the experi- 
mental curve from Ref. llSl . 



IV. ELECTRONIC STRUCTURE 
CALCULATIONS 

A. Local Density Approximation 

Since there exists no previous report on the electronic 
structure of KTi(S04)2-H20, we begin by analyzing the 
results from a non-magentic LDA calculation. In a sim- 
plified, fully ionic model each Ti'^+ ion is surrounded by 
a slightly distorted octahedron of O^^ ions. The pseudo 
octahedral coordination dictates a set of local axes for 
the conventional eg and t2g orbitals. The local coordi- 
nate system is chosen as z 1 1 c, and x and y axes are 
rotated by 45 around c with respect to the original a 
and b axes. The non-magnetic total and orbital resolved 
density of states (DOS) are collected in Fig. [2] The pre- 
sented part of the valence band is predominantly com- 
prised of Ti 3d and O 2p states belonging to TiOe oc- 
tahedra. The states belonging to sulfur (not shown) lie 
below -8 eV and are therefore well separated from the 
TiOe states. The weight close to the Fermi level (Eif) is 
mainly from the Ti t2g states which contain two electrons 
(one for each Ti in the unit cell) and are separated by a 



ligand field splitting of about 2 eV from the higher lying 
(empty) Cg states. For an octahedral arrangement of oxy- 
gen anions around a 3d transition metal cation, a 2 eV 
ligand field split is rather typical. For KTi(S04)2-H20, 
the band width of the t2g band complex is only 0.65 eV, 
about one third of the value for TiOCl (-2 eV)^. This 
difference arises from the fact that in TiOCl, the basic 
octahedral structural units TiCl204 are arranged such 
that they are corner-sharing in the a-direction and edge- 
sharing along the ^-direction, leading to a larger inter- 
action between the octahedral units and hence a larger 
t2g band width. In contrast, the TiOg octahedral units 
in KTi(S04)2-H2 are neither corner- nor edge-sharing 
and hence the smaller t2g band width. 

The degeneracy between the i2g orbitals is lifted due 
to the monoclinic symmetry of the crystal structure as 
seen in the non- magnetic band structure (Fig. [3]). There 
are 2 Ti atoms per formula unit and therefore 6 i2g bands 
close to Ep. The bands are highly dispersive along F-Y, 
X-M and XZ-MZ directions, which are along the crystal- 
lographic y-axis and remain rather flat along the other 
high symmetry directions. This implies that the main 
interaction between the Ti^+ ions is along the "double 
chain" , while sizably smaller interactions are expected 
between the adjacent double chains. The band belong- 
ing to the dxz orbital is lower in energy as compared to 
the dyz and d^z (nearly empty and larger band width) 
orbitals (also see Fig. [5]). The mixing between the Ti 
3d and the O 2p states close to Ep is less than 10% and 
similar to other systems where Ti occurs in the d^ config- 
uration. This scenario is different from the cuprates (d^) 
where 30% of the contribution to the states at Ep stems 
from O 2p. This fundamental difference in the strength 
of the hybridization between the transition metal ions 
and the oxygen ligands comes from the relative energies 
(A) of the oxygen p and the transition metal d bands. 
In cuprates, the highest (half- filled) dx^-y-2 orbital and 
the uppermost (filled) oxygen p-orbitals are rather close 
in energy (A ~ 2 eV) resulting in a strong pd hybridiza- 
tion. In titanates, on the contrary, the t2g orbitals lie 
much higher in energy than the uppermost (filled) oxygen 
p-orbitals (A > 3 eV) and therefore exhibit significantly 
less pd hybridization. Upon hole doping, the holes would 
formally appear in the oxygen p-orbitals for cuprates and 
in one of the i2g orbitals for titanates, characterizing 
them as charge-transfer and Mott-Hubbard systems, re- 
spectively. 

Experimentally, KTi(S04)2-Il20 is an insulator, but a 
metallic solution is obtained within LDA. Such metallic 
results from LDA are well known and understood to arise 
from the inadequate treatment of the strong Coulomb 
correlation of the 3(i orbitals. Therefore, the orbital de- 
pendence of the Coulomb and exchange interactions are 
taken into account in a mean field like approximation 
using the LSDA-|-[/ approach (Section IV-C). As men- 
tioned previously in Sec. II, the distorted octahedra in 
TiOCl split the i2g states into a lower lying singlet [dxy] 
and a higher energy doublet {dxz,yz)- Presuming no fur- 
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FIG. 3. (Color Online) Band structures with band character 
close to the Fermi level obtained within LDA. The bands are 
dispersive along F-Y and X-M direction - the crystallographic 
y direction along which the "double chains" propagate. There 
are two Ti atoms in the unit cell and therefore 6 t2g bands 
crossing the Fermi energy. 



ther symmetry breaking, adding correlations, the choice 
of the orbital for occupying the one unpaired electron 
of the Ti'^^ ion in TiOCl is rather straightforward: the 
singlet dxy On the other hand, for KTi(S04)2-H20 the 
pseudo octahedral ligand field splits the t2g states and re- 
moves the three fold degeneracy. Nevertheless, due to the 
relatively small energy differences, there is no clear can- 
didate for the choice of the half-filled orbital. The band 
center of the dxz orbital is around Ep and is the lowest 
lying band of the i2g complex. The center of gravity of 
dyz and the comparably broad dxy band are only slightly 
higher in energy than the dxz at 0.15 eV and 0.25 eV, re- 
spectively. Since the three t2g orbitals are quite close in 
energy, a subtle balance between the orbitals is expected 
and therefore one should carefully consider the possibil- 
ity of an orbitally ordered solution. At this juncture, it is 
rather unclear which correlated orbital will be half filled 
and thereby determine the magnetic model. Even the 
undoped, low dimensional S=l/2 cuprates that possess 
an extensive literature, where the magnetic model is gen- 
erally understood to be governed by the half- filled dx^-y^ 
orbital, can sometimes show surprises. For example, the 
CuOg octahedral environment in the insulating S=l/2 
quasi ID system CuSb2 06 is less distorted than usual, 
so that the cubic degeneracy for the eg ligand-field states 
are only slightly lifted. The energy difference between the 



,,2 and d-t 



related band centers is about 0.3 eV 



only, compared with about 2 eV for standard cuprates. 
Inclusion of correlations, changes the order and results 
in the unpaired electron occupying the d^z^-r^ orbital 
instead of the standard dx^-y^ in CuSb206.— 

We take into account the strong electronic correlations 
of the Ti 3(i states via two possible ways: (a) map- 
ping the results from LDA first to a tight-binding model 
(TBM), which in turn is mapped onto a Hubbard model. 





ti (eV) 


t2 (eV) 


Ji^PM (K) 


J^^^ (K) 


a 


dyz 
dxy 


0.047 
0.037 
0.081 


0.029 
0.028 
0.135 


31 
19 
92 


11.8 
11 
256 


0.38 
0.57 

2.78 



TABLE I. Hopping parameters (in meV) and the correspond- 
ing exchange constants (in K) from an effective one-band 
TBM. The hopping paths are indicated in the lower panel 
of Fig. [T] The transfer integral ii corresponds to hopping 
between the single chanis, while t2 corresponds to hopping 
within a single chain. A C/cff of 3.3 eV has been used to cal- 
culate the J's. The frustration ratio a = J2/J1 is collected in 
the last column. 



and subsequently to a Heisenberg model. At half fill- 
ing, when U is much larger than the bandwidth w, the 
spin degrees of freedom are well described by an S=l/2 
AFM Heisenberg Hamiltonian with an exchange interac- 
tion J^^^ ~ 4i2/C/off. (b) performing LSDA+C/d total 
energy calculations for various collinear spin configura- 
tions and mapping the energy differences on to a classi- 
cal Heisenberg model to obtain the total exchange J. At 
this juncture, the difference between the two parameters 
used to incorporate the effects of strong correlations, Ucs 
and Ud must be clarified. The former is applied to LDA 
bands which include the effects of hybridization between 
the metal atoms and ligands, while the latter is applied to 
purely atomic d-orbitals and therefore necessitates using 
different values for these two parameters. 



B. Tight Binding Model 

Though LDA fails in reproducing the insulating ground 
state of KTi(S04)2-H20 system, it still provides valuable 
information about the orbitals involved in the low en- 
ergy physics, as well as their corresponding interactions 
strengths. As mentioned previously, the strong Coulomb 
correlations favor full polarization and a detailed analy- 
sis is necessary to identify the ground state for the fully 
polarized d-orbital. Therefore, to obtain a microscopic 
picture of the magnetic interactions, we have first derived 
an effective one-band TBM for each set of the t2g bands. 
For each set of the different "active" t2g orbitals, we fit 
the corresponding LDA bands to a separate TBM. In 
second quantization formalism, the tight-binding Hamil- 
tonian can be expressed as, 

HtB = Yl *y(^IaCj<T + c]aCm) + Yl ^^^t^,, (1) 

{i,j},a i,cr 

where cj^, Cja are the usual creation and annihilation op- 
erators; a denotes the spin polarization; ti are the con- 
stant energy shifts with respect to the Fermi level; tij 
are the transfer integrals and (j,j) are the site indices. 
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FIG. 4. (Color Online) The total exchange constants J's (top 
panel) and the frustration ratio a (bottom panel) as a function 
of Ud for dxz and dy^ orbitals. The NN total exchange J\ (full 
symbols) and NNN total exchange J2 (empty symbols) do not 
vary much for the considered range of Ud (2.5 eV to 4.5 eV) 
values. The experimental results (Ref. Iisl ) are indicated by 
dashed lines. 



The magnitudes of the leading hopping integrals (the 
paths are indicated in Fig. [T]) are collected in Table H] 
The NN (between single chains) and NNN (along the 
single chains) hopping integrals ti and t2 are of the same 
order of magnitude for the dxz and dyz orbitals. On the 
contrary, the TBM fit to the broad dxy orbital leads to a 
much larger t2- All other i's beyond NNN are less than 
0.1 meV for all the three orbitals and can be neglected 
for the chain physics. The two main hopping terms are 
thus confined to interactions between the Ti sites within 
each S — 1/2 "double chain", consistent with the ex- 
perimental observations^^ of displaying low dimensional 
magnetic properties. The individual exchange constants 
are then calculated using J^^^ ^ 4t2/C/eff. For TiOCl, 
a UcS ~ 3.3 eV was shown to provide good agreement 
between calculated exchange constants and susceptibil- 
ity measurements^^. The same value of Ucs = 3.3 eV 
has been used here for KTi(S04)2-H20, the J's obtained 
are collected in Table. H] The J's for the lower energy d^z 
band and as well as the dyz band are of similar magnitude 
as compared to the experimental report. To the contrary, 
both the calculated NN and NNN J's for the dxy band 
are much larger in energy compared to the experimental 
report. This large difference in the energy scale of the 
magnetic exchange for dxy band implies that this orbital 
might be a rather unlikely choice for full polarization (a 
more clear picture emerges in the following section when 
performing LSDA-t-C/ calculations). 

Since the constructed Heisenberg Hamiltonian takes 
into account only the antiferromagnetic component of 
the couplings, inclusion of ferromagnetic (FM) compo- 
nent could change the picture described above, especially 




FIG. 5. (Color Online) Wannier functions of the three t2g orbitals of the Ti'^+ ion. Left panel: dx 
right panel: dyz- The blue and red lobes of the WF's refer to the positive and negative iso-surfaces, the circles refer to the 
atoms. The TiOe octahedra is highligted by the light grey Ti-O bonds. The perspective view is also aided by highlighting the 
crystallograhpic a ~ b plane using dashed lines. 



considering the absolute values for the AFM exchanges 
are quite small. In the proceeding section we report on 
the LSDA+t/d calculations to obtain the total exchange 
constants J = J^^ + J^""^. 



C. Density functional theory + U 

Besides obtaining estimates for all the AFM couplings 
in the system, the TBM also allows for approximating 
the number of Ti-Ti neighbors that needs to be consid- 
ered when performing the more involved and time con- 
suming LSDA+Ud supercell calculations. Since the AFM 
exchanges obtained from the TBM beyond the NNN are 
less than 0.1 meV, and because FM interactions beyond 
second neighbors should also be small, we constructed 
two supercells to obtain the values of the short-ranged ex- 
changes Ji and J2. Using different initial density matri- 
ces for the Ti 3d orbitals, one can correlate (fill the spin- 
up band with one electron and leave the spin-down band 
empty) the bands belonging to different irreducible repre- 
sentations. For KTi(S04)2-H20 we tried to spin-polarize 
each of the three t2g bands. We considered Ud values 
ranging from 2.5 - 4.5 eV.'^'^ First let us consider the dxz 
and dyz orbitals, which gave similar AFM exchange con- 
stants in our TBM. In all of our LSDA+Ud calculations, 
the scenario in which the dxz band was spin-polarized 
had the lowest energy. Spin-polarizing dyz required an 
additional energy of 350 meV per Ti^+ ion. This en- 
ergy scale is comparable to the band width of these or- 
bitals. Incidentally, all our attempts to spin-polarize d^y 
remained unsuccessful, though not unexpected when con- 
sidering our TBM results (i.e. the much too large energy 



scale of J2 which is inconsistent with experiments). We 
were also unable to stabilize different orbitally-ordered 
scenarios {i.e. one Ti ion with a spin-polarized dxz or- 
bital and the NN Ti ion with a spin-polarized dyz or- 
bital). This alludes to the fact that the (local) mag- 
netic ground state in KTi(S04)2-H20 is very likely de- 
termined by the Ti 3dxz orbital. The next question to an- 
swer is whether the exchange constants obtained for the 
dxz orbital are consistent with the experimental findings. 
We obtain effective exchange constants by performing 
LSDA+Ud calculations of differently ordered spin con- 
figurations (FM, AFM and ferrimagnetic) and map the 
energies to a Heisenberg model. Among the considered 
spin configurations, the AFM spin configuration was al- 
ways more favorable (lower total energy) than the FM. 
The exchange constants and the frustration ratio a are 
collected in Fig. |4l For comparison, we have displayed 
the values for both dxz and dyz orbitals. For the range of 
Ud values considered here, Ji"^ is comparable to exper- 
imental findings while J2^^ is larger by almost an order 
of magnitude. Comparing the total J's in Fig.|4]with the 
jAFM obtained from TBM in Table. H we can infer that 
there is a significant FM component to the first neigh- 
bor {J™) while the FM componet to the NNN (J™) is 
quite negligible. Though the J's do not vary very much 
for Ud = 2.5 to 4.5 eV, an appropriate Ud value needs 
to be chosen for comparison with experiments. Spin- 
and orbital-unrestricted Hartree-Fock calculation of the 
on-site Coulomb interaction for various transition-metal 
oxides, recommend a Ud value of 4 eV for Ti^"*" ions.— 
Using that value of Ud as a benchmark, we obtain, J.^"^ 



12 K, J° 



13.4 K and 



0.94 (Ji '^ = 9.46 K, 



J^^P = 2.8 K, and a = 0.29). The calculated value of 
c(xz = Jf ^/Jf^ is larger than the experimental value by 
a factor of 4 for the d^z orbital. The calculated ayz (= 
0.60) for the energetically unfavorable dy^ orbital is also 
somewhat larger than the experimental value. Albeit the 
S=l/2 frustrated chain magnetism in KTi(S04)2-H20 is 
established in both LDA and LSDA+Ud calculations, our 
results for a are not consistent with the experiments. 
Nonetheless, both our calculation and experiment sug- 
gest an a in the highly interesting region ( 0.2411 < a < 
1.8) of the spin-1/2 frustrated chain phase diagram. 

In the following section, we attempt to understand 
the discrepancy between our calculations and the exper- 
iment. 



V. EFFECT OF CRYSTAL WATER 

A. Wannier Functions 

The octahedron surrounding Ti'^+ consists of 5 0^~ 
ions shared with S04~ groups, and one water molecule 
(see Sec. II for a more detailed discussion). From a naive 
consideration of average equatorial and apical Ti-0 dis- 
tances, the magnetically active orbital would be assigned 
as dxy In contrast, both LDA and LSDA+Ud results 
point to dxz as the magnetically active orbital. In order 
to shed some more light on this discrepancy, we have plot- 
ted the Wannier functions for the Ti t2g complex (Fig.[5|). 
Wannier functions are essentially a real-space picture of 
localized orbitals and can be used to enhance the un- 
derstanding of bonding properties via an analysis of fac- 
tors such as their shape and symmetry. Before analyzing 
the WF's, one should keep in mind that the two near- 
est neighbor (NN) Ti atoms do not belong to the same 
chain but to the pair-chain displaced along the c-axis (see 
lower panel of Fig. [1]). By fitting (using exact diagonal- 
ization) the low temperature magnetic susceptibility, the 
recent experimental report ^^ suggests that the AF-NN 
interaction Ji is larger than the AF-NNN interaction J2 . 
Therefore, one expects to observe large tails at oxygen 
sites from the WF's bending towards the NN Ti atoms, 
as this facilitates the Ti-0-O-Ti superexchange. The Ti 
dxy WF is composed of contributions from the d^y or- 
bital as well as tails on the oxygen sites in the xy-plane, 
although these do not point towards the NN Ti atom. 
The dyz WF, on the otehr hand, has tails on all 6 oxy- 
gen sites of the TiOe octahedra, and all point towards 
the NN Ti atom. Interestingly, the d^z WF not only has 
tails on the oxygen sites bending towards the NN Ti, but 
also has tails on one of the hydrogen site belonging to the 
crystal water molecule. Such an effect is arising from the 
hybridization effect of the O and H orbitals. The effects 
of hybridization involving the H-atom in the dyz WF is 
comparatively less than the d^z, since no extended tails 
are observed on the H site. We therefore infer that the 
orientation of the water molecule plays an important role 
in deciding the magnetically active orbital for obtaining 



the ground state in KTi(S04)2-H20. More detailed anal- 
ysis follows in the next section. 



B. Rotating the water molecule 

It is normally assumed that the presence of crystal 
water in a compound, leads to only modest changes in 
the crystal field of the magnetic ion. However, it might 
not always be the case. Recent work'^^ on CuCl2 and 
CuCl2-2Il20 shows that the former is a quasi ID-chain 
compound while the latter is a classical 3D AFM, al- 
though the magnetic Cu-Cl units in both compounds are 
structurally similar. Motivated by such reports and the 
results from our WF analysis, we explore the effects of 
the orientation of the H2O molecule in KTi(S04)2-H20. 
Crystallographic data that we have used until now orient 
the H2O molecule in the xz plane. From the analysis of 
the WF's we infer that the magnetically active orbital 
(d^z) is selected by the hydrogen bonding of the water 
molecule in the xz-plane with the Ti dxz orbital. We now 
attempt to rotate the H2O molecule to the yz plane and 
observe the effects on the electronic structure. The sym- 
metry of the crystal was reduced to accommodate the 
new orientation, though the experimental O-H distance 
(0.8747 A) as weh as the H-O-H bond angle (93.33°) were 
left untouched. The resulting non-magnetic LDA band 
structure is displayed in Fig. [B] The general shape of the 
bands as well as the various k - point degeneracies are 
similar for both the xz and yz oriented H2O molecule. 
The important changes for the fictitious yz oriented sys- 
tem are: (a) the d^z band is no longer having the lowest 
band energy and is getting pushed upwards and therefore 
less occupied, (b) the dyz band is pulled downwards and is 
getting more occupied. Ergo, the orientation of the H2O 
molecule has a profound effect on the electronic structure 
and plays an important role in selecting the magnetically 
active orbital. In contrast with the LDA result, the ad- 
dition of correlations results in dxz as the lowest energy 

band, as was the case for H2O xz. That said, the 

energy difference between singly occupied dxz and dyz is 
only 35 meV, 10 times smaller than the original calcu- 
lation (cf. Sec. IV-C) in which the H2O molecule was 
oriented along the xz-plane. Calculating the J's for the 
fictitious system using the effective TBM resulted in ayz 
— 0.36 and a much larger axz = 1-11. These results 
reaffirm the the importance of the water molecule in de- 
ciding not only the magnetically active orbital, but also 
the accordingly derived exchanges and frustration ratio 
in KTi(S04)-H20. 



C. Effect of the O-H bond length 

Having observed the sensitivity of the electronic struc- 
ture to the orientation of the H2O molecule, our next item 
of interest was to check the importance of the O-H bond 
length on the electronic structure and magnetic proper- 




FIG. 6. (Color Online) Comparison of the non-magnetic band structure for KTi(S04)2-H20 using the reported a:2;-plane 
oriented (left) and the fictitious yz-plane oriented (right) H2O molecule. The insets show a small segment of the double chains 
with the two respective orientations of the H2O molecule. 



ties. Recent reports on spin-1/2 kagome lattice systems^ 
show a dramatic impact of the 0-H bond length on the 
exchange. It is well known that obtaining the correct 
0-H bond length via X-ray diffraction in a system con- 
taining heavy atoms is difficult. OH^ impurity in oxides 
have been shown to have an equilibrium bond length of 
« 1 A,^ For KTi(S04)2-H20, the reported 0-H bond 
length is 10% smaller and only 0.874 A. One possible 
reason for the larger frustration ratio a in LSDA-t-t/^ cal- 
culations as compared to experiments might arise from 
the possibly underestimated 0-H bond lengthi^ We have 
therefore allowed the 0-H bond length to relax. Keeping 
the TiOg octahedra rigid, we relaxed the H position with 
respect to the total energy and obtained an optimized 
0-H bond length of about 1 A, in accordance with the 
empirical expectations.— 1^ Recalculating the exchange 
constants using the optimized 0-H distance, we obtain 
for Ud = 4.0 eV, Jf^ = 10 K, Jf^ = 14.2 K. The frus- 
tration ratio axz — 1.4 is even larger than the previ- 
ously calculated value using the experimental 0-H bond 
length. Thus far, all of our calculations have resulted 
in exchanges that have the correct sign and compara- 
ble magnitudes to the previous experimental report, but 
nonetheless always yielded frustration ratios far greater 
than previously reported, ^^. In the following section, we 
make a final attempt to clarify the discrepancy between 
the experiment and our calculations. 



VI. TMRG CALCULATIONS 

As demonstrated in Ref. Iis!, a frustrated Heisenberg 
chain model with a = 0.29 can reproduce the experimen- 
tal magnetic susceptibility curve of KTi(S04)2-H20. The 
small a = 0.29 implies that Ji is large, while J2 is small. 
This is at odds with our LSDA-t-t/d calculations, where 
the antiferromagnetic exchange between NNN Ti atoms 
appears to be more efficient than the NN exchange. Thus, 
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FIG. 7. (Color Online) Fits to the magnetic susceptibility. 
The experimental data are adopted from Ref. [isl . Magnetic 
susceptibility of frustrated Heisenberg chains with a = 0.29 
and a — 1.5 (a = J2/J1) was simulated using TMRG. The 
simulated curves were fitted to the experiment by varying 
the fitting parameters Ji , g and the temperature- independent 
contribution xo- Inset: difference curves for both solutions. 
Ax is the difference between the simulated and the experi- 
mental value. 



in the microscopic model, a should be larger than unity. 
The question is then, whether the large a regime con- 
forms to the experimental behavior. 

The small energy scale of the leading couplings leads 
to sizable error bars for the Ji and J2 values estimated 
from LSDA+Ud calculations. To refine the model pa- 
rameters, we use the analytical expressions for the high- 
temperature part of the magnetic susceptibility of a frus- 
trated Heisenberg chain, the high-temperature series ex- 
pansion (HTSE).'^^ Typical for a local optimization pro- 
cedure, the results are dependent on the initial values. If 
we start from the Ji > J2 limit, HTSE yields Ji ~ 9.6 
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FIG. 8. (Color Online) Ground state magnetization of frus- 
trated Heisenberg chains with a = 0.29 and a = J2/J1 = 1.5 
(a = J2/J1), simulated using exact diagonalization on finite 
lattices (rings) of N = 24 spins. Note the characteristic up- 
ward bending of the a = 0.29 curve. Inset: finite-size depen- 
dence of the ground state magnetization. 



K and J2 — 2.9 K, very close to the a — 0.29 reported 
in Ref. [l5|. In contrast, if we proceed from the J2 > Ji 
regime, we obtain Ji ~ 5.4 K and J2 — 8.1 K (a = 1.5), 
in accord with our LSDA+Ud calculations. Thus, HTSE 
yields two ambiguous solutions. 

HTSE typically diverges at temperatures compara- 
ble with the magnetic energy scale (T > J). To ver- 
ify, whether both solutions agree with the experimental 
x{T) at lower temperatures, we simulate the tempera- 
ture dependence of reduced magnetic susceptibility x* 
using TMRG, and fit the resulting x*(T'/max{ Ji, J2}) 
dependencies to the experimental curve. In this way, we 
again find that besides the previously reported a — 0.29 
solution, the a = 1.5 curve with Ji = 5.4K, J2 = 8.1K, 
5 = 1.74, and Xo = 5.9-10"^ emu/mol also yields an 
excellent fit to the experimental magnetic susceptibil- 
ity (Fig. [T]). The difference curves evidence that both 
a — 0.29 and a = 1.5 provide equally good description of 
the experimental x{T) data (Fig. [71 inset). 

The coexistence of the two solutions actually mani- 
fests the inner symmetry of the frustrated chain model. 
A trivial example: the uniform chain limit can be de- 
scribed with a = -|-0 ( Ji 7^ 0, J2 = 0) as well as a = 00 
(Ji = 0, J2 7^ 0). The a = 0.29 and a = 1.5 solutions 
are also related, although in a less trivial way. To pin- 
point this relation, we revisit the phase diagram of the 
frustrated chain model. The a = limit corresponds 
to the exactly solvable gapless Heisenberg chain model. 
This GS is robust against small frustrating J2 , up to the 
quantum critical point ac — 0.2411, where a spin gap 
opens. -"^ For larger a values, the spin gap rapidly in- 



creases and reaches its maximum value A ~ 0.43 Ji at 
a ~ 0.6. Further enhancement of a reduces the spin gap. 
In the large a limit (J2 » Ji), the spin gap exhibits an 
exponential decay4i 

Thus, for a certain value of the spin gap A, there 
are two possible a values: (i) with a dominant Ji, i.e. 
from the a = 0.2411-0.6 range, and (ii) with a sizable J2 
(a = 0.6-00). Since A plays a decisive role for the shape 
of x(T), both solutions yield similar macroscopic mag- 
netic behavior. This explains the seemingly unusual fact 
that the experimental data for KTi(S04)2-H20 can be 
fitted by both a = 0.29 and a = 1.5. 

Despite the similar spin gaps, the solutions a = 0.29 
and a = 1.5 are physically different (unlike, e.g. a = 
and a = 00, that describe the same physics): spiral cor- 
relations are present in the latter case, only."^! Moreover, 
the two solutions feature substantially different correla- 
tion lengths4i Thus, the two solutions can be distin- 
guished by measuring a characteristic experimental fea- 
ture ( "smoking gun" ) . 

For spin systems, a measurement of magnetization 
isotherms is technically simple, but very informative, es- 
pecially for systems with weak magnetic couplings. Since 
the magnetic field linearly couples to the S^ compo- 
nent of the spin, a magnetization curve reflects the en- 
ergy of the lowest lying state in each S^ sector. This 
often suffices to distinguish between ambiguous solu- 
tions. For instance, HTSE for the Ji — J2 square lat- 
tice system BaCdVO(P04)2 yielded, besides the frus- 
trated solution with an AFM J2, also a nonfrustrated 
solution with AFM J2. However, the frustrated sce- 
nario was clearly underpinned by M{H) measurements 4^ 
In a recent work, M{H) measurements for A2CUP2O7 
[A — Li,Na) resolved previous controversies concerning 
the magnetic dimensionality of these compounds4^ 

We argue that for the frustrated Heisenberg chain 
model, the characteristic behavior of magnetization on 
the verge of saturation can be used to distinguish between 
different scenarios. In particular, the a = 0.29 magneti- 
zation curve exhibits a well pronounced upward bending, 
while only a feeble bending is visible in the a = 1.5 GS 
magnetization (Fig. |S]). Another relevant quantity is the 
saturation field: 
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where S'^ax corresponds to the fully polarized state. The 
energies are estimated using exact diagonalization for fi- 
nite chains of A^ = 32 spins. Adopting Ji, J2, and g val- 
ues from the HTSE fits, we obtain TJgat — 16.4 T and 
_ffsat — 18.7T for a = 0.29 and a = 1.5, respectively. Both 
values of saturation field lie in the experimentally acces- 
sible field range. A somewhat problematic point could 
be the low energy scale of KTi(S04)2-H20, which ren- 
ders the typical measurement temperature of ~1.5K as 
relatively high, hence the states with different S^ could 
be substantially mixed. Still, the a = 0.29 magnetiza- 
tion isotherms will retain fingerprints of the characteris- 
tic bending. Therefore, we believe that a high-field (up 



to '^20 T) measurement of a magnetization isotherm will 
be an instructive and decisive experiment to distinguish 
between the a = 0.29 and a = 1.5 scenarios. 



VII. SUMMARY 

In conclusion, we have studied the electronic structure 
of KTi(S04)2-H20 in detail using DFT based calcula- 
tions. The resuhs of both the TBM and LSBA+Ud calcu- 
lations confirm beyond doubt the low dimensional nature 
of the material with NN and NNN exchanges Ji and J2 
confined to the double-chains running along the 6-axis. 
We also confirm the AFM nature of the exchanges, con- 
sistent with the experimental report, with the Ti-3dxz 
orbital being the magnetically active one, holding the sin- 
gle unpaired electron of the Ti^+ ion. The magnitude of 
the calculated J's are of the right order compared to the 
experiment, though we observe a strong dependence to 
the t2g orbital choice. Notwithstanding the small energy 
scale (w 10 K) of the system, we are able to obtain ac- 
curate results from our DFT calculations. Additionally, 
we observe a surprisingly large spin-lattice coupling with 
respect to the crystal water related degrees of freedom, 
which in turn manifests itself by playing an important 



role in determining the ground state of the system. This 
feature is clearly elucidated by calculating the Wannier 
functions, which show the effects of hydrogen bonding 
to the corresponding t2g orbital which is oriented in the 
same plane as the crystal water molecule. Using the ex- 
perimental position for hydrogen, we obtain a frustration 
ratio a w 1.0 and a value of a ~ 1.4 upon relaxing the 
hydrogen position in the crystal lattice. Both these val- 
ues are larger than the experimental value acxp = 0.29. 
In order to understand the origin of this discrepancy be- 
tween the experiment and our calculations, we simulated 
the temperature dependence of the susceptibility using 
both the small and large values of a. Due to an intrinsic 
symmetry of the Ji — J2 frustrated chain model, we show 
that both values of a provide equally good fits to the 
experimental curve. Consequently, we calculated magne- 
tization curves as a means to unambiguously distinguish 
the two solutions and show two features which can be 
used to identify the appropriate a that defines the mag- 
netic ground state of KTi(S04)2-H20. Hence, we sug- 
gest performing high-field magnetization measurements 
on this system and as well as a.c. susceptibility exper- 
iments (to obtain the size of the spin-gap) at very low 
temperatures. 
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